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Computation  of  Discrete  Hole  Film 
Cooling  Flow  Using  the  Navier-Stokes  Equations 
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ABSTRACT 

An  analysis  and  computational  procedure  are  described  here  for  predicting 
flow  and  heat  transfer  which  results  from  coolant  injection  through  a  single 
row  of  round  holes  oriented  at  an  angle  to  a  flat  surface  with  the  injection 
and  free  stream  velocity  vectors  coplanar.  The  present  method  solves  the  compres¬ 
sible  Navier-Stokes  equations  and  utilizes  "zone  embedding",  surface-oriented 
coordinates,  interactive  boundary  conditions,  and  an  efficient  split  LBI  scheme. 
The  approach  treats  the  near-hole  flow  region  where  the  film  cooling  flow  is 
initially  established.  A  sample  laminar  flow  calculation  is  presented  for  an 
injection  angle  of  45  degrees  and  a  ratio  of  normal  injection  to  free  stream 
velocity  of  0.1.  Selected  results  are  compared  with  previous  calculations  for 
the  normal  injection  case.  Although  the  present  results  do  not  include  heat 
transfer  predictions,  details  of  the  interaction  between  injectant  and  main 
stream  flow  near  the  hole  exit  are  in  qualitative  agreement  with  experimental 
observations  for  other  flow  conditions. 
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INTRODUCTION 


To  achieve  higher  turbine  efficiencies,  designers  have  been  forced  to  deal 
with  high  turbine  inlet  temperatures  which  may  exceed  structural  limitations  of 
available  materials  used  for  blading  unless  some  method  of  cooling  these  blades 
is  introduced.  Various  cooling  schemes  have  been  devised  to  ensure  that  the 
structural  integrity  of  the  blades  remain  intact  when  exposed  to  a  high  tempera¬ 
ture  environment  for  prolonged  periods  of  time.  One  of  the  more  promising 
schemes  is  discrete  hole  film  cooling,  wherein  cooling  air  is  injected  through 
either  a  single  row  or  multiple  staggered  rows  of  holes,  to  provide  a  protective 
layer  of  cool  air  between  the  blade  surface  and  hot  mainstream  gases.  Discrete 
hole  film  cooling  is  better  suited  to  the  high  temperatures  of  advanced  high 
performance  engines  than  the  current  convert ively  cooled  blades  found  in  commer¬ 
cial  engines.  A  large  number  of  parameters  can  be  varied  to  obtain  a  configura¬ 
tion  which  optimizes  heat  transfer  and  aerodynamic  characteristics,  including 
hole  shapes  and  patterns,  injection  angles  and  rates,  surface  curvature,  coolant 
temperatures,  and  mainstream  boundary  layer  characteristics.  Since  experimental 
determination  of  optimal  configurations  is  costly,  a  computational  procedure 
which  could  be  used  to  screen  alternative  conf igurations  without  experimental 
testing  would  be  of  considerable  value.  An  analysis  and  computational  procedure 
is  described  here  for  predicting  flow  and  heat  transfer  which  results  from 
coolant  injection  through  a  single  row  of  round  holes  oriented  at  a  non-normal 
angle  to  a  flat  surface.  The  work  discussed  herein  is  an  extension  of  a  previous 
effort  which  considered  only  normal  injection  through  a  single  row  of  round  holes. 
This  study  provides  an  improved  understanding  of  the  film  cooling  process, 
particularly  in  the  near-hole  region,  and  represents  the  second  step  toward 
development  of  a  computer  program  capable  of  predicting  heat  transfer  and 
aerodynamic  loss  levels  associated  with  discrete  hole  film  cooling  on  gas 
turbine  blades. 
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PREVIOUS  WORK 


Previous  experimental  work  on  film  cooling  flow  has  consisted  mainly  of 
flow  visualization  studies  and  measurements  of  film  cooling  effectiveness  down¬ 
stream  of  the  coolant  injection  hole  for  a  variety  of  flow  conditions.  For 
example,  Cooladay  and  Russell  [1]  performed  a  flow  visualization  study  of  dis¬ 
crete  hole  film  cooling  in  three  different  hole  configurations,  to  obtain  a 
better  understanding  of  flow  behavior  associated  with  coolant  injection.  They 
considered  injection  from  discrete  holes  in  a  three-row  staggered  array  with 
five-diameter  spacing  and  three-hole  angles:  (1)  normal  to  the  surface, 

(2)  slanted  30  degrees  to  the  surface  but  aligned  (coplanar)  with  the  free 
stream,  and  (3)  slanted  30  degrees  to  the  surface  and  45  degrees  laterally  to 
the  free  stream.  For  flow  conditions  typical  of  gas  turbine  applications, 

Colladay  and  Russell  observed  that  normal  injection  is  subject  to  flow  separa¬ 
tion  even  at  low  injection  rates,  that  slanted  injection  works  well  except  at 
high  injection  rates,  and  that  compound  si  ant ed-yawed  injection  is  less  susepti- 
ble  to  separation  hut  tends  to  promote  entrainment  of  free  stream  fluid  (which 
augments  heat  transfer)  and  to  increase  secondary  flow  losses.  These  experimental 
observations  emphasize  the  importance  of  interaction  between  the  injected  fluid 
and  the  free  stream. 

Kadotani  and  Goldstein  [2]  have  performed  experimental  measurements  for 
discrete  hole  film  cooling  through  a  row  of  holes  having  three-diameter  lateral 
spacing,  inclined  at  35  degrees  to  the  injection  surface,  and  coplanar  with  the 
free  stream.  They  considered  a  variety  of  flow  conditions  and  varied  boundary 
layer  thickness  and  injection  rates,  as  well  as  Reynolds  number  and  turbulence 
properties.  They  measured  lateral  and  streamwise  distributions  of  film  cooling 
effectiveness  downstream  of  injection,  and  also  measured  mean  velocity  and 
temperature  in  a  transverse  plane  normal  to  the  free  stream  and  located  at  the 
downstream  edge  of  the  injection  holes.  Kadotani  and  Goldstein  deduced  from 
their  measurements  that  reversed  flow  is  present  near  the  injection  hole  even 
for  injection  rates  as  low  as  0.35,  and  that  the  size  and  strength  of  the 
reversed  flow  increases  with  increasing  boundary  layer  thickness.  They  also 
found  that  the  approaching  boundary  layer  thickness  significantly  affects  the 
flow  behavior  and  lateral  variation  of  film  cooling  effectiveness,  particu¬ 
larly  near  the  injection  hole. 
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THE  PRESENT  APPROACH 


The  present  study  focuses  on  the  near-hole  flow  region  where  the  film 
cooling  flow  is  initially  established  and  where  local  hot  spots  may  occur. 

Since  the  flow  field  surrounding  injection  may  involve  separation  and  reversed 
flow  relative  to  the  free  stream,  and  since  directional  and/or  other  properties 
of  the  flow  which  might  permit  simplifying  assumptions  are  not  clearly  evident, 
the  present  approach  is  based  on  numerical  solution  of  the  three-dimensional 
compressible  Navier-Stokes  equations.  The  governing  equations  in  Cartesian 
coordinates  are  cast  into  strong  conservation  form  using  a  Jacobian  transforma¬ 
tion  to  obtain  the  governing  equations  in  a  general  nonor thogonal  coordinate 
system.  The  coordinate  system  used  in  the  present  study  is  a  skewed  elliptical- 
polar  coordinate  system  which  is  not  aligned  with  the  free  stream  flow.  However, 
all  solid  surfaces  and  computational  boundaries  lie  along  coordinate  lines  in 
the  transformed  coordinate  system. 

In  selecting  the  computai tonal  domain,  a  "zone  embedding"  approach  is 
adopted  whereby  attention  is  focused  on  a  subregion  of  the  overall  flow  field 
in  the  immediate  vicinity  of  the  discrete  hole  coolant  injection.  At  the 
curved  boundaries  located  within  the  free  stream  region,  interactive  boundary 
conditions  which  permit  inflow  and  outflow  of  shear  layers  and  the  compressible 
inviscid  free  stream  are  derived  from  an  assumed  flow  structure.  Since  details 
of  the  coolant  velocity  distribution  at  the  hole  exit  are  not  known  a  priori 
and  presumably  will  depend  on  an  interaction  between  the  coolant  and  the  free 
stream  flow,  the  computational  domain  is  chosen  to  include  the  flow  region 
within  the  coolant  hole  as  well  as  that  exterior  to  the  hole.  Including  the 
flow  region  within  the  coolant  hole  permits  the  interaction  between  the  coolant 
and  free  stream  flows  and  its  important  influence  near  the  hole  to  be  determined 
as  part  of  the  final  solution,  without  the  need  for  simplifying  assumptions. 

The  governing  equations  are  solved  by  an  efficient  and  noniterative  time- 
dependent  "linearized  block  implicit"  (LBI)  scheme. 
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METHOD  OF  SOLUTION 

Zone  Embedding  and  Interactive  Boundary  Conditions 

The  skewed  elliptical  polar  coordinate  system  (r,6,z')  and  the  computational 
grid  used  in  the  present  calculations  are  shown  in  Fig.  1,  where  z'  is  the  sheared 
coordinate  variable.  The  polar  coordinates  (r,6)  in  the  surface  plane  (z/R=0) 
vary  from  circular  at  the  origin  to  elliptical  at  the  edge  of  the  coolant  hole 
to  circular  again  at  the  outer  boundary  of  the  domain  (r/R=5).  The  latter 
choice  is  not  a  necessity,  and  the  outer  transverse  boundary  locaLion  (at  0  =  tt/2) 
may  be  selected  to  determine  the  coolant  hole  spacing.  The  injection  angle 
chosen  for  the  present  case  is  45  degrees.  Above  the  hole  the  coordinate 
system  skewing  angle  was  gradually  changed  from  45  to  90  degrees  at  the  outer 
boundary  (z/R=5)  allowing  symmetry  conditions  to  be  easily  applied,  so  that  the 
flow  represented  is  that  of  opposing  coolant  holes  located  on  parallel  flat 
plates  with  spacing  h/R=10.  Velocity  boundary  conditions  at  no-slip  and  sym¬ 
metry  surfaces  are  straightforward  and  self-explanatory.  The  remaining  condi¬ 
tion  at  the  no-slip  surface  in  the  plane  of  the  hole  (z/R=0)  was  taken  as  the 
normal  (z-cartesian)  momentum  equation,  and  within  the  coolant  hole  9p/3n  =  0 
was  applied,  where  n  is  tiie  normal  direction  at  any  surface.  The  latter  condi¬ 
tion  at  a  no-slip  surface  is  correct  to  order  Re  1  for  viscous  flow  at  high 
Reynolds  number. 

The  treatment  of  inflow  and  outflow  conditions  is  an  obstacle  to  be  overcome 
within  the  zone  embedding  approach.  This  problem  was  successfully  treated  for 
the  normal  injection  case  by  Kreskovskv,  Briley  and  McDonald  [3]  and  thus  their 
general  approach  is  followed  in  the  present  study.  At  the  computational  boundaries 
located  within  the  free  stream  region,  interactive  inflow-outflow  boundary 
conditions  tire  derived  from  an  assumed  flow  structure  and  physical  approxi¬ 
mations  following  Ref.  3.  The  assumed  inflow  structure  consists  of  a  uniform 
potential  flow  core  velocity  U^,  two-dimensional  estimates  of  the  boundary  layer 
displacement  thickness  s*(x)  on  the  flat  surface,  and  finally  an  estimate  of  the 
blockage  correction  factor  B(x).  The  boundary  layer  displacement  thickness 
6*(x)  is  estimated  from  the  known  development  of  flow  between  parallel  flat 
plates  (e.g.  Kreskovsky  and  Shamroth  16]  and  references).  Given  this  estimate 
of  6*(x),  the  boundary  layer  thickness  is  approximated  by  the  Blasius  relation¬ 
ship  d ( x )  =  5/1.72  6*(x)  and  the  blockage  factor  is  given  by  B(x)  =  [H/2R~6*(x)]  ^ 

Finally,  boundary  layer  velocity  profile  shapes  f  (z/6),  0-f-l  are  defined  from 
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von-Karman  Pohlhausen  polynomial  profiles  and  the  boundary  layer  thickness  d(x). 
The  initial  velocity  vector  U  at  t=0  is  defined  by 

U(r,0,z)  =  UjB  (x)  f  [z/S{x)j  U) 

At  the  inflow  boundary,  a  “two-layer"  boundary  condition  is  devised  such 

that  stagnation  pressure  P  is  fixed  at  the  free  stream  reference  value  (P  ) 

o  or 

in  the  core  flow  region  (z>6),  and  the  cartesian  velocity  u  is  set  by  u  = 
ue(x,t)  f(z/6)  for  z  -  <5.  The  local  free  stream  velocity  ug  is  chosen  consis¬ 
tent  with  the  free  stream  total  pressure  Pq,  total  temperature  T  ,  and  the  local 
static  pressure  pg  at  a  specified  point  on  the  boundary,  such  that  the  initial 

distribution  of  u  is  maintained: 
e 

ue(x,n  =  ue(x,1  =  0)-  f(p0,  T0,  pe)  (2) 

and  the  function  f  follows  from  isentropic  relations.  The  static  pressure  at 

the  point  in  question  is  determined  as  part  of  the  solution  and  is  updated 

after  each  time  step.  The  latter  application  of  the  two-layer  boundary  condition 

is  a  variation  of  the  condition  used  by  Ref.  3,  which  was  required  because  of 

the  skewed  nature  of  the  computational  domain.  The  remaining  inflow  conditions 

2  2  2  2 

are  v  =  0,  9  v/9n  =  0,  and  3  c  /9n  g(6,z)  where  g  is  the  value  of  this  quantity 

p2-  - 

at  t  =  0  with  0^  defined  as  (1-B  which  is  its  value  from  the  potential 

flow  corrected  for  estimated  blockage. 

For  outflow  conditions,  second  derivatives  of  each  computed  velocity 
component  are  set  to  zero  and  the  static  pressure  is  imposed  and  updated  after 
each  time  step  from  an  assumed  form  which  varies  linearly  with  x  from  a  fixed 
value  at  the  downstream  limit  of  the  flow  region  (r/R=5,  P=0)  to  a  value  at 
6  =  Jl/2  which  is  consistent  with  the  upstream  total  pressure  and  the  local  free¬ 
st  ream  dynamic  head.  At  the  inflow  boundary  within  the  coolant  hole,  the 
normal  velocity  component  is  given  a  specified  velocity  profile  distribution 

with  uniform  core  region,  and  remaining  velocity  components  are  extrapolated 

2  2 

using  zero  second  derivative.  The  remaining  condition  is  9  p/9n  =  0. 

The  foregoing  interactive  inflow-outflow  boundary  conditions  are  designed 
to  permit  the  mass  flux  through  the  computational  domain  to  adjust  to  both  the 
imposed  downstream  static  pressure  and  to  viscous  losses  present  in  the  flow, 
while  maintaining  a  specified  flow  structure  based  on  physical  assumptions 
consistent  with  the  flow  problem  under  consideration.  Various  refinements  in 
the  interactive  boundary  conditions  are  possible,  such  as  including  the  effect 
of  local  pressure  gradient  on  boundary  layer  growth  and  profile  shape. 
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Dif ferencing  Procedures; 

The  differencing  procedures  used  are  a  straight  forward  adaptation  of  those 
used  by  Briley  and  McDonald  [7]  in  cartesian  coordinates  for  flow  in  a  straight 
duct.  A1 though  the  present  procedure  can  ultimately  be  used  to  compute  local 
heat  transfer  distributions,  the  present  study  has  focused  on  methodology  for 
treating  the  fluid  flow  near  and  within  an  injection  hole.  Accordingly,  for 
the  present,  the  fluid  stagnation  enthalpy  is  assumed  constant  for  reasons  of 
economy.  The  definition  of  stagnation  enthalpy  and  t  lie  equation  of  state  for  a 
perfect  gas  can  then  he  used  to  eliminate  pressure  and  temperature  as  dependent 
variables,  and  solution  of  the  energy  equation  is  unnecessary.  The  continuity 
and  three  cartesian  momentum  equations  are  solved  in  the  (r,P,z')  coordinate 
system  with  density  and  the  cartesian  velocity  components  (u,v,w)  as  dependent 
variables.  Three-point  central  differences  were  used  for  spatial  derivatives, 
and  second-order  artificial  dissipation  terms  are  added  as  in  [7]  to  prevent 
spatial  oscillations  at  high  cell  Reynolds  number.  This  treatment  lowers  the 
formal  accuracy  to  first  order  hut  does  not  seriously  degrade  accuracy  in 
representing  viscous  terms  in  thin  shear  layers.  Analytical  coordinate  trans¬ 
formations  were  used  to  provide  increased  resolution  of  critical  portions  of  the 
flow  field.  These  coordinate  t  ransforKi.it  ions  serve  to  concentrate  grid  points 
near  the  flat  surface,  near  the  hole  surface,  near  the  hole  exit,  and  in  the 
region  just  downstream  of  the  hole.  Derivatives  of  geometric  data  were  deter¬ 
mined  numerically  for  use  in  the  strong  conservation  form  of  the  difference 
equa t  i ons . 

Spl  it  I.R1  Al  gor  i  t  hm 

The  numerical  algorithm  used  in  the  consistently-split  "linearized  block 
implicit"  (LR1)  scheme  developed  hv  Briley  and  McDonald  (7,8]  for  systematic 
use  in  solving,  svsiems  of  nonlinear  pa rahol i c -hyperbol i c  partial  differential 
equations  (DDK's) .  To  illustrate  the  algorithm,  let 

{(j)0*  '  -  <£n)/At  =  /3D(c£n  +  ')  +  (I  ~/3)  D(cjbn)  (3) 

appro:  imat  e  a  system  of  t  i  rr.f  -d<  pendent  nonlinear  PDE '  s  (centered  about  tn+f  At) 
for  the  vector  <}  of  dependent  variables,  were  D  is  a  mult  idimensional  vector 
spatial  differential  operator,  and  t  is  a  discretized  time  variable  such  that 
At  =  t ^  —  t n .  A  local  time1  linearization  (Taylor  expansion  about  <n)  is  intro¬ 
duced,  and  this  selves  to  define  a  linear  differential  operator  L  such  that 
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D  (<£n*')  =  D  (<jtn )  +  Ln(<£n+I-  £")  +  0(  At2) 


(4) 


Eq.  (3)  can  thus  be  written  as  the  linear  system 

(I  -£A11_n)(<£n+,-<£n)  =  AID (<£n) 


(5) 


The  multidimensional  operator  L  is  divided  into  three  "one-dimensional"  sub¬ 
operators  L  =  +  L2  +  (associated  here  with  the  three  coordinate  directions), 

and  Eq.  (5)  is  split  as  in  the  scalar  development  of  Douglas  and  Gunn  [9]  and  is 
written  as 


(I  -  /3  At  L  1°)  ( (£*  -  c £n)  =  AtD(c£") 

(6a) 

(I  -/9AtL2n)(<£**  -  <£n)  =  <£«-</>" 

(6b) 

(1  -J3 AtL3n)(<£***  -  </>"): 

(6c) 

4>n+'-- $***  +  0 ( At 3 ) 

(6d) 

If  spatial  derivatives  appearing  in  L  are  replaced  by  three-point  difference 
formulas,  then  each  step  in  Eqs.  (6a-c)  can  be  solved  by  a  block-tridiagonal 
"inversion".  Eliminating  the  intermediate  steps  in  Eqs.  (6a-d)  results  in 

(I  -  A  t  L  ,n )  (  3  -y3A1L?n)(I  -|SAtL3n)ic/>n4l-^n)  -  A1  D  (<*>")  (7) 

3 

which  approximates  Eq .  (5)  to  order  At  .  Complete  derivations  are  given  in 
[7,8].  It  is  noted  that  Beam  A  Warming  [10]  have  reformulated  this  algorithm 
as  a  widely-used  "delta  form"  approximate  factorization  scheme  whose  two-level 
form  is  identical  to  Eq.  (6). 

Finally,  it  is  noted  that  a  special  treatment  for  radial  boundary  conditions 
is  needed  at  the  polar  axis  r  =  0  of  the  coordinate  system  being  used.  As  a  result 
of  the  splitting  of  Eq.  (3),  a  separate  (i.e.,  "multivalued")  radial  boundary 
condition  is  needed  for  each  r-implicit  step  of  the  split  algorithm.  As  dis¬ 
cussed  in  Ref.  3  symmetry  conditions  are  valid  at  r=0  and  provide  a  usable 
radial  boundary  condition  for  the  implicit  step  which  corresponds  to  6=11/2. 
Consequently,  for  each  z  value,  the  r-implicit  steps  are  reordered  so  as  to 
solve  the  6  =  11/2  step  first,  using  symmetry  conditions.  This  defines  implicit 
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values  for  the  dependent  variables  at  r=0  which  in  turn  are  used  to  specify 
function  boundary  conditions  for  the  remaining  r-implicit  steps  (O^Jl/2),  whose 
order  is  immaterial. 

Computed  Results 

To  illustrate  the  application  of  the  present  approach  to  the  discrete  hole 
film  cooling  problem,  a  sample  calculation  has  been  computed  for  laminar  flow 
with  coolant  injected  at  45  degrees  to  the  free  stream.  Although  the  present 
results  are  limited  to  laminar  flow  with  the  coolant  injection  velocity  coplanar 
with  the  free  stream,  this  approach  may  be  extended  to  treat  turbulent  flow, 
heat  transfer,  and  a  coolant  injection  velocity  not  coplanar  with  the  free  stream. 
In  the  flow  considered,  the  ratio  of  average  coolant  velocity  to  average  main 
stream  velocity  is  0.1,  the  free  stream  Reynolds  number  based  on  duct  width  is 
400,  and  the  reference  Mach  number  is  0.3.  The  Reynolds  number  based  on  hole 
diameter  is  80.  The  thickness  6/R  of  the  approaching  boundary  layer  is  1.8. 

The  solution  converged  in  about  50  time  steps  and  with  a  14  x  19  x  15  grid 
required  about  24  minutes  of  CDC  7600  CP  time. 

Selected  results  from  this  sample  calculation  are  shown  in  Figs.  2-7. 

For  comparison  purposes  certain  results  from  the  normal  injection  case  pre¬ 
sented  in  Ref.  3  are  reproduced  in  Figs.  8-11,  A  vector  plot  of  velocity  in 
the  vertical  symmetry  plane  (6=0,11)  passing  through  the  center  of  the  coolant 
hole  is  shown  in  Fig.  2.  It  can  be  seen  that  most  of  the  coolant  flow  leaves 
the  hole  near  its  downstream  edge,  and  for  this  particular  set  of  flow  condi¬ 
tions,  the  coolant  flow  does  not  penetrate  beyond  the  mainstream  boundary  layer. 
Comparing  the  velocity  vectors  with  those  from  the  normal  injection  case. 

Fig.  8,  it  can  be  seen  that  the  penetration  of  the  jet  into  the  main  flow 
is  less  in  the  45°  injection  case,  as  expected.  Contours  of  velocity  components 
comprising  the  velocity  vectors  in  Fig.  2  are  shown  in  Fig.  3.  It  can  be  seen 
that  the  streamwise  component  is  not  greatly  disturbed  by  the  coolant  injection 
at  this  ratio  of  injectant  to  free  stream  velocity  (0.1).  However,  it  is  clear 
from  the  contours  of  vertical  velocity  that  there  is  significant  local  inter¬ 
action  between  the  main  stream  and  coolant  stream  near  the  hole  exit,  and  that 
this  interaction  affects  the  coolant  flow  within  the  hole.  It  can  be  seen  from 
these  contours  that  the  coolant  velocity  profile  becomes  highly  distorted  at  the 
hole  exit  plane  (z/R=0).  The  streamwise  and  vertical  velocity  contours  from 
the  normal  injection  case  are  shown  in  Fig.  9,  from  which  it  is  evident  the 
primary  difference  in  the  two  cases  occurs  in  the  vicinity  of  the  upstream  edge 
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of  the  coolant  hole.  There  in  the  45°  injection  case  the  vertical  component  of 
velocity  varies  rapidly  near  the  upstream  side  of  the  hole  in  a  manner  similar 
to  that  occurring  along  the  downstream  side  as  the  hole  exit  plane  is  approached. 
This  behavior  was  not  present  in  the  normal  injection  case.  Fig.  9b,  and  pre¬ 
sumably  the  phenomenon  predicted  in  the  45°  injection  case  is  due  to  the  occur¬ 
rence  of  a  low  pressure  region  near  the  upstream  edge  of  the  coolant  hole  at  the 
exit  plane.  Further  verification  of  this  behavior  would  require  recalculation 
with  a  refined  computational  mesh.  The  present  approach  provides  a  means  for 
detailed  prediction  of  this  local  interaction,  by  including  a  portion  of  the 
hole  within  the  computed  flow  region.  It  is  also  noted  that  the  present 
results  did  not  predict  a  reversed  flow  region  downstream  of  the  hole,  probably 
due  to  the  relatively  low  injection  rate  and  the  injection  angle. 

A  further  indication  of  the  coolant  injection  flow  is  given  in  Fig.  4, 
where  contours  of  the  vertical  velocity  (normal  to  the  flat  surface)  are  shown 
for  selected  horizontal  planes.  Note  that  the  location  of  the  coolant  hole 
at  z/R=0  is  shifted  to  the  left  in  the  figures,  since  the  grid  is  skewed  in  the 
z-direction.  A  vector  plot  of  velocity  in  a  horizontal  plane  one  grid  point 
above  the  hole  exit  (z/R=0.15)  is  shown  in  Fig.  5,  and  this  figure  shows  the 
extent  of  lateral  flow  deflection  present  in  the  computed  results.  The  cor¬ 
responding  vector  plot  from  the  normal  injection  case  is  shown  in  Fig.  10, 
where  comparison  reveals  the  expected  reduced  lateral  deflection  of  the  main 
stream  flow  in  the  45°  injection  case.  Contours  of  the  streamwise  (u)  and  trans¬ 
verse  (v)  components  of  velocity  shown  in  Fig.  6  at  z/R=0.15  further  illustrate 
the  influence  of  the  coolant  injection  on  the  mean  flow  boundary  layer.  Finally, 
vector  plots  of  secondary  velocity  in  vertical  planes  perpendicular  to  the  main 
flow  direction  are  shown  in  Fig.  7  for  several  streamwise  locations  upstream, 
within,  and  downstream  of  the  hole.  These  plots  provide  an  indication  of  the 
secondary  flow  and  fluid  entrainment  induced  by  the  coolant  injection.  Again, 
this  behavior  is  qualitatively  similar  to  the  secondary  flow  patterns  observed 
in  the  normal  injection  case.  Fig.  11. 

A  problem  which  makes  a  detailed  assessment  of  the  predictions  difficult  is 
the  occurrence  of  pressure  oscillations  in  the  computations  near  the  edge 
of  the  coolant  hole  exit  plane.  Although  the  solution  remains  stable  in  their 
presence  the  oscillations  do  make  a  quantitative  assessment  of  the  results  in 
the  vicinity  of  the  coolant  hole  difficult.  The  oscillations  appear  to  emanate 
from  the  edge  of  the  hole  which  is  a  reentrant  corner  in  the  computat ional 
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(r,  z')  plane  for  each  value  of  6.  These  corner  points  are  presently  treated  in 
the  split  LBI  framework  using  the  "double-value"  concept,  whereby  values  for  each 
variable  at  the  corner  are  associated  with  the  direction  from  which  the  corner  is 
approached.  In  the  present  case  the  velocities  are  identically  zero  at  the  cor¬ 
ner,  so  only  density  and  pressure  are  multivalued.  This  approach  tends  to 
reduce  the  region  of  solution  contamination  around  the  corner,  however,  an 
alternate  approach  to  the  double  values  may  be  necessary  to  reduce  the  oscillations 
and  obtain  quantitatively  accurate  results. 

Summary  and  Concluding  Remarks 

A  computational  procedure  has  been  described  for  predicting  flow  and  heat 
transfer  which  results  from  coolant  injection  through  round  holes  oriented  at 
45  degrees  to  a  flat  surface.  The  procedure  solves  the  compressible  Navier- 
Stokes  equations  without  simplifying  assumptions,  and  includes  portions  of  the 
flow  field  both  within  the  coolant  hole  and  exterior  to  the  hole.  Zone  embed¬ 
ding  and  iteractive  boundary  conditions  are  used  to  achieve  economy  by  limiting 
the  solution  domain  to  the  near-hole  flow  region,  where  simplifying  assumptions 
are  undesirable  because  of  local  iteraction  and  reversed  flow.  A  sample  coarse- 
mesh  calculation  has  been  obtained  for  laminar  adiabatic  flow  with  a  ratio  of 
normal  injection  to  free  stream  velocity  of  0.1.  This  calculation  demonstrates 
the  general  capability  of  the  present  procedure  for  treating  coolant  injectant 
flows  without  making  simplifying  assumptions  in  the  near-hole  flow  region. 

The  results  obtained  for  the  45  degree  injection  case  are  qualitatively 
reasonable,  and  the  differences  observed  between  this  case  and  the  normal 
injection  case  appear  to  be  physically  plausible.  A  more  detailed  evaluation 
of  the  procedure  and  its  capabilities  would  require  fine-mesh  calculations  and 
consideration  of  nonadiabatic  flows  and  other  coolant  velocity  ratios.  The 
present  procedure  can  also  be  extended  to  treat  compound  injection  through 
holes  not  coplanar  with  the  free  stream,  by  modifying  the  treatment  of  the 
coordinate  system  singular  points  (r=0)  and  by  solving  the  equations  over  the 
full  angular  range  (0-0-211)  with  periodic  boundary  conditions. 
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(a)  St  Teamwise  Velocity  Component 


(b)  Vertical  Velocity  Component 


Fig.  3  -  Contours  of  Velocity  in  W-rtica}  Symmetry  Plane 
Bisecting  Coolant  Hole,  45°  Injection. 
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(a)  Streamwise  Velocity  Component 
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Fig.  6  -  Contours  of  Velocity  in  Hoi  i.'ontal  pl-ne  at  x/K  0.15, 
45°  1 n  j  ec t  ion. 
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Secondary  Velocity  Vectors  in  Vertical  Planes  Perpendicular  to  Main 
Flow  Direction,  45°  Injection. 
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Fie.  9  -  Contours.  of  Velocity  in  Vertical  fyr.etiv  Plane  Bisecting 
Coolant  I'ole.  N'cir.a.l  Inject  icn. 


